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ABSTRACT: 

The  US  Army  Research  Laboratory  (ARL)  is  currently  conducts  tests  on  anti-ballistic 
armor  for  thilitary  uses.  This  research  is  concerned  with  determining  the  limit  velocity  (vL)  of 
different  target  penetrator  combinations.  The  limit  velocity  is  the  highest  velocity  a  penetrator 
can  have  without  penetrating  the  targe.  Unfortunately,  penetration  processes  are  highly  complex 
and  an  effective  first  principles  derivation  of  vL  has  not  beat  discovered.  Estimation  of  vL  is 
therefore  done  empirically.  Furthermore,  ballistics  tests  can  be  very  expensive,  resulting  in  a 
small  size  sample  with  which  to  perform  statistical  data  analysis. 

There  are  two  ballistics  testing  methods  commonly  used  to  estimate  vL.  The  Jonas- 
Lambert  method  involves  measuring  the  residual  velocity  of  the  projectile  after  perforation.  The 
bisection  method  or  V30  simply  evaluates  the  perforation  without  residual  velocity.  The  second 
method  is  significantly  less  expensive. 

Simulation  is  used  to  model  both  of  the  common  ballistics  testing  methods  as  well  as 
several  new  approaches  to  ballistics  testing.  The  results  are  evaluated  and  compared  for 
statistical  significance  and  accuracy.  This  work  suggests  that  the  bisection  method  is  more 
accurate  when  sample  size  is  small.  This  discovery  could  provide  considerable  cost  savings  to 
ballistics  testing  at  ARL. 
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INTRODUCTION 


Problem 

The  US  Army  Research  Laboratory,  ARL,  is  currently  responsible  for  testing  the  armor 
the  military  uses.  This  research  is  concerned  with  determining  the  limit  velocity,  VL,  for 
different  types  of  armor.  The  limit  velocity  is  the  fastest  velocity  a  round  being  shot’ at  a  certain 
type  of  armor  can  have  without  penetrating  that  armor.  Knowledge  of  the  limit  velocity  will 
provide  a  metric  for  the  optimality  of  a  given  type  of  armor.  In  addition,  the  limit  velocity  can 
be  used  to  establish  the  maximum  velocity  an  enemy  projectile  can  have  without  allowing  a 
penetration. 

Unfortunately,  there  have  been  no  effective  first  principle  derivations  for  limit  velocity 
developed,  therefore  empirical  testing  is  used  to  find  an  accurate  prediction.  However,  this 
experimental  testing  is  extremely  expensive  due  to  the  number  of  real  rounds  fired,  the  armor 
that  is  used  for  the  testing,  and  any  equipment  used  for  measuring  data.  Therefore,  a  method 
must  be  developed  to  minimize  the  sample  size  of  shots  fired  while  still  allowing  a  true 
prediction  for  limit  velocity.  •  *• 

Currently  there  are  two  techniques  for  this  testing:  The  Jonas  -  Lambert  Method  which 
t  11868  thelspeed  of  the  round  both  before  and  after  penetration  and  Bisection  Method,  orVx,  which 
does  not.  The  latter  is  much  less  expensive.  This  study  will  compare  the  currentmefeods  and 
conclude  which  one  is  optimal  in  terms  of  accuracy,  cost,  and  statistical  significance  . 
Additionally,  different  techniques  that  are  not  being  used  will  be  developed.  These  new  methods 
wiH  alst>  be  compared  with  fee  two  current  techniques  to  determine  if  a  “better”  system  for 
testing  exists.  ^  .. 

Background  Information  ;  .  „ 

The  limit  velocity  is  defined  as  the  highest  striking  velocity  a  piece  of  armor  can 
withstand  without  allowing  a  complete  penetration  by  a  ballistic  round.  Thus,  each  piece  of 
armor  will  have  a  different  limit  velocity  based  upon  the  type  of  round  fired  at  it  This 
relationship  is  seen  in  Figure  1  below: 

v„| 


Figure  1:  Vsv  VK 


Vs 


VR  repr  esents  fee  residual  velocity,  the  speed  of  the  round  after  penetration  occurs,  and  Vs 
is  the  striking  velocity,  fee  speed  of  the  round  before  contact  with  the  armor.  Every  striking 
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velocity  that  has  a  positive  residual  velocity  represents  a  complete  penetration.  As  the  residual 
velocity  becomes  smaller,  the  striking  velocity  approaches  the  limit  velocity.  The  limit  velocity 
is  the  fastest  striking  velocity  that  has  a  residual  velocity  of  zero. 

Finding  the  limit  velocity  in  a  perfect  world  with  the  above  information  would  be  fairly 
simple.  All  that  would  be  required  is  a  function  that  describes  the  relationship  between  the 
striking  velocity  and  residual  velocity.  The  root  of  this  function  could  be  solved  for  and  this 
would  be  the  limit  velocity.  The  impact  dynamics  in  this  type  of  testing,  however,  do  not  act  in  a 
perfect  way.  An  extensive  search  of  the  literature  reveals  that  an  accurate  relationship  between 
the  striking  velocity  and  the  residual  velocity  has  not  yet  been  derived.  Figure  2  shows  a  scatter 
plot  of  a  typical  striking  velocity  versus  residual  velocity  relationship.  In  the  real  world  there  is  a 
significant  amount  of  error  involved  and  there  exists  a  zone  of  mixed  results. 


The  circled  region  in  Figure  2  is  referred  to  as  a  “zone  of  mixed  results.”  As  a  result  of 
unknown  factors,  some  test  shots  fired  at  lower  striking  Velocities  can  produce  larger  residual 
velocities.  This  makes  finding  a  relationship  between  the  residual  and  striking  velocities  very 
difficult  The  armor  is  not  flawless  and  thus  some  striking  velocities  that  are  above  the  limit 
velocity  will  not  completely  penetrate  the  armor.  Likewise,  some  striking  velocities  below  the 
limit  velocity  will  completely  penetrate  the  armor.  This  zone  is  real  and  will  not  disappear. 
Every  test  researched  contained  this  error  to  some  extent.  For  the  purposes  of  this  research,  this 
zone  is  assumed  to  be  small  enough  to  not  greatly  affect  the  results  and  will  not  be  included. 

METHODS 

Bisection,  or  Vx,  Method 

There  are  two  common  methods  for  testing  limit  velocity.  Neither  technique  has  been 
proven  “optimal”  yet.  The  cheaper  of  the  two  methods  is  The  Bisection,  or  V„ ,  Method.  Based 
on  the  Intermediate  Value  Theorem,  this  technique  looks  for  the  limit  velocity  by  treating  the 
Vs  v  Vr  plot  as  a  continuous  and  differentiable  function,/  The  function  is  defined  by  an  interval 
[a,b]  with  jXa)  mdffb)  of  opposite  signs.  However,  for  this  curve,  zero  values  are  considered  to 
have  negative  values  because  there  are  no  negative  residual  velocity  values.  By  the  Intermediate 
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midpoint  value.  If  a  “  mZt ZlS^S?  ‘ S,Wved  “d  ‘  “*  ah°'  is  hdcen  at  this  tower 
2£  oer,  ,„,  is 


Figuw  3 :  The  Bisection  Method 

zets^seest^ 

His  method  is  rehuively  cheap  bec^f  “  **”*»*■ 

only  coneemed  with  whether  or  tS^S^  c^T^^  ™»  ““  “ 

and  the  pieces  of  amor  used  cover  the^riotfthe  entire  ^  '  P"“  “oh  ** 

The  Jonas  -  Lambert  Method 

““ti,  ss  “ 

concerned  with  whether  a  complete  penetration  ocSX  ThC  bC1”8 

—*  A  graphic  4SZ  StSSh  J©  S'  *  ^  Wi“  ^  ^  Bmi' 
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Vs 

Figure  4:  The  Jonas-Lambert  Curve 


It  can  be  seen  in  Figure  4  above  that  the  root  of  the  function  is  equal  to  the  limit  velocity. 
The  equation  of  this  line  is  of  the  form: 

a(V/-Vl'/'  +  e,Vl<Vs 

0,0  £VS  £Vl 

This  equation  is  derived  from  The  Law  of  Conservation  of  Energy.  The.  derivation  and 
an  explanation  of  it  is  included  in  Appendix  A. 

This  method  requires  at  least  three  shots  where  the  striking  velocity  is  greater  than  the 
limit  velocity  to  effectively  estimate  the  limit  velpcjty,  Three  shots  are  required  in  order  to  give 
three  separate  equations  which  in  turn  allow  for  solving  the  three  unknown  parameters, 
excluding  the  error  term.  After  defining  The  Jonas  -  Lambert  Curve,  the  limit  velocity  can  be 
found  by  simply  finding  its  root. . „ .  . .  . 

An  efficient  application  of  The  Jonas  -  Lambert  Method  involves  taking  three  shots  that 
completely  penetrate.  These  three  shots  will  provide  foe  necessary  d«r»  to  solve  for  all  the 
unknown  parameters  except  the  error  term.  The  calculated  limited  velocity  is  then  used  as  foe 
next  shot’s  striking  velocity.  If  a  complete  penetration  occurs  on  this  next  shot,  foe  variables  are 
re-solved  and  the  next  shot  will  be  at  foe  latest  calculated  value  of  foe  limit  velocity..  If  a 
complete  penetration  does  not  occur,  re-shoot  at  a  higher  velocity  until  a  penetration  occurs. 

Once  a  penetrations  occurs,  continue  with  the  method.  Shots  will  continue  to  be  fired  until  an 
answer  of  desired  accuracy  is  found. 


Vs  v  Vg  Relationships 

The  first  method  is  looking  for  a  relationship,  similar  to  the  one  used  in  foe  Jonas  - 
Lambert  Method,  between  the  residual  and  striking  velocities.  If  an  accurate  function  which 
describes  the  residual  velocity  based  on  the  striking  velocity  can  be  found,  the  root  of  this 
function  will  provide  a  more  accurate  estimate  of  the  limit  velocity. 

This  method  requires  foe  same  testing  techniques  and  data  and  thus  will  be  gimilar  in 
price  to  The  Jonas  -  Lambert  Method.  The  difference  is  that  common  relationships  are  used  to 
describe  the  data  rather  than  some  complex  equation.  This  method  looks  at  linear,  exponential, 
power  and  logarithmic  relationships. 
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Golden  Ratio  Method 


The  Golden  Ratio  Method  is  an  extension  of  The  Bisection  Method  and  is  very  similar. 
This  method  uses  The  Golden  Ratio,  a  ratio  that  appears  consistently  throughout  nature,  as  the 
way  of  determining  the  bracket  used  rather  than  simply  using  a  halving  technique.  Like  the 
Bisection  Method,  this  method  is  only  concerned  with  whether  or  not  the  round  penetrates  the 
piece  of  armor  or  not  Therefore,  the  price  of  this  method  only  deals  with  how  many  shots  are 
taken.  Additionally,  the  same  initial  bracket  of  +/-  50  m/s  will  be  used  to  remain  consistent  with 
the  evaluation  of  The  Bisection  Method. 

The  algorithm  for  this  method  starts  with  an  initial  bracket  of +/-  50  m/s  and  shoots  at  the 
middle.  The  initial  bracket  is  divided  into  an  upper  and  a  lower  61 .8  %  (the  golden  ratio) 
portion.  If  the  first  shot  penetrates,  the  lower  portion  is  used  as  the  second  bracket  and  if  a 
penetration  does  not  occur,  the  upper  portion  is  used.  The  second  shot  then  shoots  at  the  middle 
of  this  second  bracket.  Once  again,  if  a  penetration  occurs,  the  lower  61.8%  of  the  second 
bracket  is  used  as  the  third  bracket,  and  if  no  penetration  occurs,  the  upper  61 .8%  portion  is  used. 
This  iterative  technique  continues  until  a  limit  velocity  of  desired  accuracy  is  determined.  This 
process  is  seen  in  Figure  5. 


h* 


h* 


Figure  5:  The  Golden  Ratio  Method 

The  lower  of  the  two  overlapping  lines  represents  the  upper  61 .8%  of  the  bracket  and  the 
higher  of  the  two  is  the  bottom  6 1 .8%.  The  first  shot,  labeled  “1 ,”  bisects  the  initial  bracket  in 
half.  The  “X”  represents  the  limit  velocity  and  thus,  the  first  shot  is  a  penetration.  The  lower 
portion  is  then  used  as  the  next  bracket.  The  “2”  marks  the  second  shot  which  halves  the  bottom 
61 .8%  of  the  bracket  This  shot  is  a  penetration  as  well  and  thus  the  lower  portion  is  used  for  the 
third  bracket.  The  second  green  line  down  becomes  the  new  bracket  and  foe  third  shot  will  then 
half  this  new  bracket.  This  iterative  method  continues  until  foe  limit  velocity  is  found  to  within  a 
desired  level  of  accuracy. 

Residual  Energy  v  The  Angle  of  the  Projectile  Relationships 

This  method  deals  with  trying  to  find  a  relationship  between  the  angle  of  foe  projectile, 
with  respect  to  foe  armor,  when  it  strikes  the  target  and  the  residual  energy.  The  residual  energy 
is  foe  energy  of  foe  projectile  after  it  penetrates  the  armor.  As  foe  round  hits  foe  armor,  energy  is 
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absorbed  by  the  armor,  the  round  breaks  up,  and  energy  is  contained  in  all  of  the  round 
fragments. 

Using  the  principles  of  conservation  of  energy,  an  equation  is  derived  that  gives  the  limit 
velocity  as  a  function  of  this  residual  energy  and  the  striking  mass  of  the  round.  The  equation  is 
of  the  form: 


The  derivation  of  this  is  in  Appendix  B.  Therefore,  if  a  relationship  can  be  found  which 
describes  this  residual  energy,  E,  based  on  the  angle,  then  the  limit  velocity  can  be  found. 

This  method  investigated  may  be  costly  because  data  must  be  gathered  on  the  striking 
angle  of  the  round.  The  particular  factors  of  the  angle  in  this  method  are  the  pitch,  yaw  and 
resultant  of  the  projectile  as  it  strikes  the  armor. 

Residual  Energy  v  Vg  Relationships 

This  method  investigates  a  relationship  between  the  residual  energy  and  the  residual 
energy.  The  goal  is  to  find  a  function  that  describe  foe  residual  energy  based  on  the  residual 
velocity.  The  root  of  this  function  will  yield  the  residual  energy  at  the  limit  velocity*  which  can 
then  be  used  to  solve  for  the  limit  velocity. 

This  method  requires  some  costly  data  collection 'as  well.  The  residual  and;atriking 
velocity  as  well  as  foe  striking  and  residual  mass  of  the  round  must  be  measured:to  define  the 
function.  The  residual  velocity  is  obviously  needed  because  it  is  the  independent  yariable.  The 
other  pieces  of  data  are  needed  to  solve  for  foe  residual  energy.  The  residual  energy  is  found 
using  foe  equation  of  foe  form: 

..  ..  -^MSVS.  MtVg. 5=2^ 

The  derivation  of  this  equation  is  located  in  Appendix  B. 

RESULTS 

Bisection  or  Vx  Method 

The  Bisection  Method  can  produce  a  limit  velocity  to  any  desired  accuracy.  The 
algorithm  increases  accuracy  with  every  iteration  performed.  To  test  this  method,  a  code  in 
Microsoft  Excel  was  developed  that  takes  successive  shots  until  foe  desired  accuracy  is 
produced.  For  foe  test,  sample  limit  velocities  were  used  from  data  given  from  ARL.  A  few 
examples  of  this  code  is  in  Appendix  C.  Table  1  shows  the  results  of  all  foe  tests  conducted. 
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trial 

shots  for  a 
tolerance  of  5  m/s 

shots  for  a 
tolerance  of  1  m/s 

1 

3 

5 
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3 

F  4 
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4 
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5 

2 

5 

6 

r  3 
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2 

6 

8 

3 

3 

9 

2 

6 

10 

3 

6 

11 

3 

6 

12 

4 

6 

average 

shots 

£916666667 

5.5 

•  Table  1:  Bisection  Results 


It  can  be  seen  in  Table  Lthatthis,method  produces  a  result  within  5  m/s  of  the  true  limit 
velocity  with  2.9  shots  on  average;  .A  point  of  interest  with  these  results  is  the  consistency; 
-There  are  no  outliers  jn  the  data.  The  mpstjshots  eyer.  needed  to  reach  this  accuracy  of  5  tn/s  is 
only 4 shots.  *  m!, 

The  Jonas  -  Lambert  Method  ... 

•  *  ;  "  r  ; 

:»■ 

,  TheJonas  -  Lambert  Method  does  not  accurately  and  completely  describe  the  behavior 
of  a  projectile  penetrating  a  piece  of  armor.  For  example,  two  shots  that  have  die  exact  same 
staking  velocity  may  result  in  different  residual  velocities.  This  is  the  main  reason  why  this 
method  has  problems.  To  test  this  method,  data  sets  obtained  from  ARL,  which  include  all  the 
information  needed  to  use  The  Jonas  -  Lambert  Method,  were  used.  This  is  in  Appendix  D. 

The  Jonas  -  Lambert  Method’s  estimate  was  compared  with  the  actual  limit  velocity  The  results 
are  shown  in  Table  2. 


complete 

penetrations 

shots 

a 

P 

vL 

estimate 

(m/s) 

Vl  act 
(m/s) 

error 

(m/s) 

4 

6 

0.89 

2.59 

1365.684 

1373 

7.316 

3 

5 

0.9 

2.49 

1082.559 

1088 

5.441 

3 

5 

0.86 

1303.539 

1324 

20.461 

3 

8 

1 

2.99 

1164 

27.875 

3 

9 

0.88 

2.59 

1178.759 

1239 

60.241 

2 

9 

0.97 

5.69 

1341.986 

1355 

wrzim 

average 

7 

22.39133 

Table  2;  Jonas  -  Lambert  Results 
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,  _  With  the  data  supplied  by  ARL,  The  Jonas  -  Lambert  Method  was  used  to  estimate  the 
limit  velocity  for  each  shot  To  arrive  at  the  limit  velocity  estimation  seen  in  Table  2,  every  shot 
from  each  data  set  is  used  to  calculate  a  limit  velocity  and  then  the  average  from  the  entire  data 
set  gives  the  final  estimate.  A  table  showing  this  data,  along  with  a  sample  algorithm  for  this 
method,  is  in  Appendix  E. 

The  data  shows  that  this  method  on  average  takes  7  shots  to  reach  a  limit  velocity  that  is 
accurate  to  within  22.4  m/s.  However,  this  data  can  be  deceiving.  First,  this  average  error,  is 
just  that,  an  average.  Table  2  shows  that  outliers  existed  where  the  method  was  only  accurate  to 
within  60  m/s.  The  average  number  of  shots  is  a  better  representation  of  the  data,  because  the 
most  shots  ever  taken  in  the  testing  was  9.  However,  there  is  no  substantial  relationship  between 
the  number  of  shots  taken  and  then  accuracy  of  the  estimate.  For  example,  5  shots  produced  an 
estimate  accurate  to  within  5.4  rri/s  and  9  shots  gave  two  different  less  accurate  estimates:  60.2 
m/s  and  13.0  m/s.  This  is  different  from  what  is  expected;  in  theoiy,  as  the  number  of  shots 
increase,  the  accuracy  should  increase  as  well. 

Vs  v  V/f  Relationships 

*  This  method  was  stfipped  before  any  results  were  calculated  because  the  findings.were  so 
inaccurate.  The  problem  is  that  the  Vs  v  VR  data  does  not  fit  into  a  “normal”  (logarithmic,  , 

exponential,  power,  or  linear)  function.  Out  of  these  functions,  the  best  relationship  was  the 
logarithmic.  The  results  for  this  method  are  shown  in  Table  3. 


trial 

shots 

penetrations 

vL 

actual 

(m/s) 

VL 

estimate 

(m/s) 

error 

(m/s) 

1 

6 

4 

1366 

1343.46 

22.54043 

2 

5 

3 

1083 

1074.171 

291.8289 

3 

5 

3 

1304 

1247.197 

118.8031 

5.333333 

144.3908 

Table  3:  V5v  V*  Results 


Only  three  sequences  of  this  method  occurred  because  the  results  were  already 
unacceptable.  As  you  can  see  the  average  error  was  already  with  144.4  ra/s  and  one  estimate 
was  off  by  almost  300  m/s.  ThiB  is  unacceptable  because  ballistidans  can  usually  guess  to 
within  100  m/s.  Sample  calculations  for  this  method  are  in  Appendix  F. 

Golden  Ratio  Method 

This  method  was  tested  in  the  same  manner  as  The  Bisection  Method.  Microsoft  Excel 
was  used  to  perform  iterations  of  the  algorithm  until  a  desired  estimate  was  calculated.  A  few 
sequences  of  the  test  is  in  Appendix  G.  The  results  of  this  method  are  shown  in  Table  4. 
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trial 

shots  for  a 
tolerance  of  5  m/s 

shots  for  a 
tolerance  of  1  m/s 

l 

2 

6 

2 

3 

5 

3 

5 

6 

4 

4 

4 

5 

3 

5 

6 

3 

3 

7 

2 

4 

8 

4 

4 

9 

3 

6 

10 

4 

4 

11 

3 

5 

12 

3 

5 

average 

shots 

3.25 

4.75 

Table  4:  Golden  Ratio  Results 


.  The  results  from  Table4  show  that  this  method  takes  3.3  shots  on  average  to  achieve  an 
estimate  within  5  m/s.  It  should  be  noticed  that,  like  the  bisection  method,  this  method  is  Very 
consistent  and  has  no  outhera.  A  point  of  interest  about  this  method  in  comparison  to  the  ^ 
Bisection  method  is  thediflerence  between  achieving  an  estimate  to  within  1  m/a  versus  5  m/s 

ItlT*  i0tS/°r/T5f  Cf5lden  111,110  Mcthod  the, Bisection  Method  to  achieve  an 

estimate  to  within  5  m/s  but  less,  shots  for  1  m/s.  .  , 


Residual  Energy  v  The  Angle  of  the  Projectile  Relationships 

^estimates  of  the  limit  velocity  are  presented  because  no  significant  relationships  coul 

7S'?'  “S  **  ofthe  PK>j«tile.  Out  of  the  three  factors,  th 
Hot?18  ^  rcIadonsiuP  but  it  was  still  very  insignificant  and  inconsistent 

The  ddta  supplied  by  ARL  along  with  the  calculated  residual  energy  is  located  in  Appendix  H 
1  contains  a  few  sample  regressions  between  beta  and  the  residual 
energy.  As  mentioned  earlier,  yaw  produced  the  most  significant  results,  therefore  only  these 
bccauaethe  other  attempts  proved  fruitless.  The  exponential  trmd  line  is 
displayed  because  it  was  the  only  fimetion  that  could  be  possibly  used  for  each  dat*  set. 


Residual  Energy  v  Vs  Relationships 


for  «s  to  find  a  relationship  between  the  residual  energy  and  the 

a,relationshiP  significance  was  found  to^tkue 

wifr  the  test  Appendix  J  shows  a  few  ofthe  sample  regressions  between  these  two  variables  as 
well  as  sample  calculations.  Table  5  shows  the  results  of  this  method. 
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shots 


i 


i 


trial 

_ 1_ 

_ 2_ 

_ 3_ 

_ 4. 

5 


complete 

penetrations 


relationship 
significance 
(r  squared) 


VL 

estimate 
(m/s) 


VL 

actual 

(m/s) 


error 

(m/s) 


0.693 


1345.095 


1373 


27.905 


0.729 


1076.384 


1088 


11.616 


0.983 


1245.647 


1324 


78.353 


0.83 


995.6095 


1164 


168.3905 


0.252 


978.7987 


1239 


260.2013 


6.6 


109.2932 


Table  5:  Residual  Energy  v  VR  Relationships  Results 


This  method  produces  very  fickle  results.  Some  tests  show  that  an  error  of  only  1 1 .6  m/s 
can  be  reached  by  only  taking  5  shots,  but  others  result  with  a  260.2  m/s  error  with  9  shots.  The 
average  results  show  this  method  can  produce  an  estimate  with  an  error  of  109.3  m/s  with  6.6 
shots.  This  test  is  automatically  discarded  because  the  error  is  larger  than  the  bracket  that  the 
ballisticians  can  initially  guess. 

*  CONCLUSION 

•  .  The  Bisection  Method  is  the  optimal  method  to  find  the  limit  velocity.  This  method  is 
..  -considered  the  best  method  based  on  cost,  accuracy,  and  reliability.  The  Bisection  Method  is  the 
cheapest  way  to  find  the  limit  velocity.  Not  only  docs  this  test  not  require  any.  expensive  data 
collection  equipment  but  it  requires  the  least  amount  of  shots  to  run  the  test  .  More  importantly, 
this  method  provided  fire  most  accurate  results.  For  a  desired  accuracy  of  5  m/«;.'The  Bisection 
Method  is  able  to  achieve  these  results  with  every  test  Lastly,  this  method  is  the  most  reliable  as 
well.  Outliers  do  not  exist  in  this  method;  plus  dr  minus,  two  shots,  this  method  always  returns 
an  accurate  answer. 

The  Golden  Ratio  Method  places  a  close  second  to  The  Bisection  Method  for  the  same 
three  reasons  discussed  above.  However,  The  Bisection  Method  is  able  to  outperform  the 
Golden  Ratio  Method  for  the  number  of  shots  to  produce  an  estimate  within  5  m/s  and  thus  The 
Bisection  Method  is  chosen  as  the  better  of  the  two. 
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recommendation 

The  following  algorithm  should  be  applied  for  finding  the  limit  velocity. 

1 .  -Estimate  the  limit  velocity  and  shoot  at  this  velocity. 

2.  -If  a  complete  penetration  occurs,  use  this  velocity  as  the  right  limit. 

-If  a  complete  penetration  does  not  occur,  use  this  velocity  as  the  left  limit. 

3  thefirstshot6^  61101  *  +  °f  ”  75  of  Ieft  or  limit  depending  on  the  result  of 


wfft  ,h“  “““ comp,rcd  *i,h  !»“»■ 

-If  the  same  result  for  the  second  shot  occurs  as  compared  with  the  first  shot,  take  the 
next  shot  at  +  or  -  25  m/s  of  the  last  shot  in  order  to  ensure  a  +  or-  lOOm/s  bracket  is 
achieved.  If  the  opposite  result  from  the  previous  shots  still  does  not  occur,  continue 
taking  shots  at  +  or  -  25  m/s  until  a  complete  penetration  or  no  penetration  occurs 
dep”ldl”8  on  what  result  is  needed.  Once  a  bracket  of  a  complete  penetration  and’  no 
penetration  is  developed,  proceed  with  the  bisection  algorithm. 


83 


WORKS  CONSULTED 


Lambert,  I.  P.  and  G.  H.  Jonas,  Towards  Standardization  in  Terminal  Ballistics  Testing:  Velocity 
Representation  (BRL  Report  No.  1852).  Aberdeen  Proving  Grounds,  MD:  USA  Ballistic 
Research  Laboratories,  1976 

Leonard,  Wendy.  The  Effect  of  Nose  Shape  on  Depleted  Uranium  (DU)  Long- 

Rod  PenetratorsfARL  —  TR  —  15 05).  Aberdeen  Proving  Grounds,  MD:  Army  Research 
Laboratory,  1997. 

McCulloh,  Ian  A.  and  William  deRosset.  Statistically  Valid  Estimation  of  Limit  Velocity  in 

Ballistics  Testing  with  Small  Sample  Size.  West  Point,  NY:  Department  ofMatbematical 
Sciences,  West  Point  and  US  Army  Research  Laboratory.  . 


84 


Appendix  A 

Hie  Jonas  -  Lambert  Derivation 

^  *»  to 

energy  of  the  romd^en  tbe  Sf^bdty^u^100^'  Therefore  E  is  ctJual  to  the  kinetic 


E-W' 


in  the  test  ***  not  change  at  all  points 

rZ-V+V 

1.  •  •  * 

.  * 

,  ,  ■*' '  “^S*  +^t*).*  +•*:. 

“  '  .  r; ,  *.$  ,  ;  .♦  .,  . 

mere  u  an  error  term  added  because  some  oftfae  ’ 

masses  are  not  exact  Jonas  —  t  *u ■■smnpiions  used  for  the  cancellation  of 

help  describe  wyowwrTte^rifam^6"  ^  2  •  top’s  and  adds  an  afytha  tom  to*  / 
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Appendix  B 


Residual  Energy  Derivation 

Start  with  the  conservation  of  energy  equation.  The  energy  of  the  striking  round,  the  left  side  of 

.*  e5ft,on* wlU  °?ual  of  the  residual  round  plus  the  energy  lost  in  the  armor  and 

with  fragments  of  the  round,  the  right  side  of  the  equation. 

2msf's  ~ +  E 

The  energy  of  the  main  piece  of  the  residual  round  is  then  subtracted  from  both  sides  There  is 
no  summatron  attached  to  the  kinetic  energy  of  the  residual  round  because  the  term  only  deals 
with  tiie  main  residual  mass  that  penetrates  the  armor.  Now  the  term  on  the  right  side  of  the 

Si  811  °f  **  magy  ‘S  l08t  a*mOT  through  *=  kinetic  energy  in 


-  “  | msV,2 =%£ 


At  the  limit  velocity,  the  residual  velocity  of  the  round  will  be  0.  The  equation  is  changod  to  this 
instance.  The  residual  velocity  goes  to  zero  and  the  striking  velocity  becomes  the  limit  velocity. 


j 


jmsKt  -£E 

»  •  •;  ••  a  Jt.  “  ;  *.  * 


This  equation  is  manipulated  to  give  the  limit  velocity  as  a  function  of  the  energy  lost  and  the 
mass  of  the  striking  round. 
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Appendix  C 


Sample  Bisection  Method  Tests 


1300 

Tolerance 

5 

VL 

actual 

1366 

rtaht  limit 

1400 

left 

ism 

Success? 

shot  1 

1350 

1300 

1400 

Fail 

shot  2 

1375 

' 

1350 

1400 

Fall 

shot  3 

1362.5 

1350 

1375 

Success 

Trial  1 


SEMI 

1000 

Tolerance 

5 

vL 

actual 

1083 

1100 

BFhwbsb 

left 

EH 

Success? 

shot  1 

1050 

fflar^T 

1100 

IZlflEMi 

shot  2 

1075 

1050 

1100 

Fall 

shot  3 

■■■RE5Z1 

. 

.  1075 

1100 

il 

Success 

Trial  2 


tmsmmm 

1300 

Tolerance 

5 

VL 

actual 

1304 

riaht  limit 

mm 

■ 

- 

■left 

Success? 

shot  1 

1350 

1300 

1400 

FaH 

shot  2 

1325 

1300 

1350 

FaB 

shot  3 

1312.5 

1300 

1325 

FaH 

shot  4 

1306.25 

1300 

1312.5 

Success 

Trial  3 


left  limit 

1100 

Tolerance 

5 

vL 

actual 

1136 

1200 

MUM 

left 

ElvUHK 

Success? 

shot  1 

1150 

1100 

1200 

FaB 

shot  2 

1125 

‘ 

1100 

1150 

FaB 

shot  3 

1137.5 

1125 

1150 

Success 

Trial  4 


note  -  all  velocities  are  measured  in  m/s. 
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Appendix  D 


Supplied  Data 


Va _ 

-Vs 

a 

P 

Vl  actual 

293 

1405 

0.89 

2.59 

1373 

360 

1363 

198 

1378 

847 

1437 

360 

1123 

0.9 

2.49 

1088 

359 

1108 

91 

1089 

318 

1330 

0.88 

2.39 

1324 

383 

1344 

473 

1371 

824 

1285 

1 

2.99 

1184 

772 

1232 

538 

1188 

618 

1254 

0.88 

1239 

419 

1288 

561 

1285 

801 

1381 

■SI 

1046 

1384 

0.87 

5.69 

1355 

655 

1370 

note  -  all  velocities  are  measured  in  m/s. 
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Appendix  E 

Jonas  -  Lambert  Calculations 


IRH 

\m 

m 

P 

Vl  estimate 

IE9M 

293 

360 

196 

647 

1405 

1383 

1378 

1437 

1392.25619 

1360.596933 

1373.381776 

1336.50248 

26.57185 

-5.08741 

7.697431 

-29.1819 

in 

360 

359 

91 

1123 

1108 

1089 

2.49 

1087.672754 

1072.183109 

1087.822589 

5.11327 

-10.3764 

5.263105 

1082.559 

1088 

318 

383 

473 

1330 

1344 

1371 

0.86 

2.39 

SI 

■  SfSOT 

1324 

824 

772 

538 

1285 

1232 

1168 

1 

2.99 

■ 

■  SSSa 

23.21946 

-15.7428 

-7.4787 

1136.125 

618 

419 

581 

801 

■reji 

■  SI 

0.88 

2.59 

1137.727821 

1249.301445 

1177.459531 

1150.54833 

-41.0315 

70.54216 

-1.29975 

-28.211 

1178(759 

J 

! 

1 

1046 

855 

1384 

1370 

0.97 

5.69 

1318.302238 

1365.669192 

-23.6835 

23.68348 

1341  {986 

i 

;  :  '1355 

,» *  t 

note  -  all  velocities  are  measured  in  m/s. 
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Appendix  F 


Vs  v  Vjt  Calculations 


am 

mm 

Vi 

actual 

shot  1 

1406 

293 

1366 

shot  2 

1383 

360 

shot3 

1378 

196 

shot  4  ' 

1437 

647 

1300  1300-  1400  1420 


Vs 


Bn 

K9H 

fam 

shot  1 

1123 

1083 

shot  2 

■KE1 

nr a 

shot  3 

MtfT  l 

91 

1000  1090  1100  1110  1120  1130 

Vs 


I  Vl  estimate  1  1074.171  1 
Trial  2 


all  velocities  are  measured  in  m/s. 
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Appendix  G 


Sample  Golden  Ratio  Method  Tests 


note  -  all  velocities  are  measured  in  m/s. 
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Appendix  H 


Data  For  The  Residual  Energy  v  The  Angle  of  the  Projectile  Relationships 


— — S—  mR - Ms - E _ flamma  aloha  beta 

66.33  1376  4.37  196  62862549  1  58  1  c . 

6626  m  «38  360  63011343  0.34  0.25  025 

“g  llll  293  66151776  1  6  1-25  f 

67236540  0,89  ns 

VL  actual 
1373 

66  29  US  926  360  41143322  0.79  0.25  0  75“ 

S2  106  769  3S0  40182488  0.58  0.25  0  5 

- - - -91  39326178  0.34  an  nos 

1088 

65.78  1330  6.55  318  57847940  0.25  0 

65.79  1344  7.56  383  58864638  0.56  0  25  Vs 

-fg - _ 473,  61274391  07  Vt  “ 

1324 

65.75  1285  10.9  824  50583603  1  1  n 

1232  10-96  772  48670417  0.56  0  5  0  25 

— ffff - ~§f - 1^2 - 536  43835053  •  1.51  i  s  I'll 

1164 

MM  1288  7.24  410  84151005  “58  f ! 

6603  1285  0.09  581  .51401021  0.75  0  75  0 

- — - in _ Sir  58615812  0  80  p” 

1239 

V  ll.v 


.<<  ' 


note  -  all  velocities  are  measured  in  m/s. 
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Appendix  I 

Sample  Regressions  For  The  Residual  Energy  v  The  Angle  of  the  Projectile  Relationships 


note  -  gamma,  alpha,  and  beta  represent  fee  resultant  angle,  pitch  and  yaw  of  fee  round, 
respectively 
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EvVr 

y  =  2615.2Ln(x)- 46416 


57000000  58000000  58000000  60000000  61000000  62000000 


E 


IDs 

E  at  Vl 
estimate 

VL  estimates 

VL  actual 

VL 

estimate 

error 

65.78 

51061786.44 

1245.994329 

1324 

1245.647 

78.35257 

65.79 

51061786.44 

65.88 

1245.048315 

Tria 

13 

note  -  all  velocities  are  measured  in  m/s. 
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